home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_syr2k.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  20.4 KB  |  828 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_syr2k (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int uplo = 121;
  14.    int trans = 111;
  15.    int N = 1;
  16.    int K = 2;
  17.    float alpha = -0.3f;
  18.    float beta = 0.1f;
  19.    float A[] = { -0.635f, 0.805f };
  20.    int lda = 2;
  21.    float B[] = { 0.773f, 0.375f };
  22.    int ldb = 2;
  23.    float C[] = { 0.616f };
  24.    int ldc = 1;
  25.    float C_expected[] = { 0.174988f };
  26.    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  27.    {
  28.      int i;
  29.      for (i = 0; i < 1; i++) {
  30.        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1622)");
  31.      }
  32.    };
  33.   };
  34.  
  35.  
  36.   {
  37.    int order = 101;
  38.    int uplo = 122;
  39.    int trans = 111;
  40.    int N = 1;
  41.    int K = 2;
  42.    float alpha = -0.3f;
  43.    float beta = 0.1f;
  44.    float A[] = { -0.396f, -0.131f };
  45.    int lda = 2;
  46.    float B[] = { -0.603f, -0.288f };
  47.    int ldb = 2;
  48.    float C[] = { -0.434f };
  49.    int ldc = 1;
  50.    float C_expected[] = { -0.20931f };
  51.    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  52.    {
  53.      int i;
  54.      for (i = 0; i < 1; i++) {
  55.        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1623)");
  56.      }
  57.    };
  58.   };
  59.  
  60.  
  61.   {
  62.    int order = 102;
  63.    int uplo = 121;
  64.    int trans = 111;
  65.    int N = 1;
  66.    int K = 2;
  67.    float alpha = -0.3f;
  68.    float beta = 0.1f;
  69.    float A[] = { -0.085f, -0.444f };
  70.    int lda = 1;
  71.    float B[] = { 0.936f, 0.752f };
  72.    int ldb = 1;
  73.    float C[] = { -0.64f };
  74.    int ldc = 1;
  75.    float C_expected[] = { 0.184069f };
  76.    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  77.    {
  78.      int i;
  79.      for (i = 0; i < 1; i++) {
  80.        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1624)");
  81.      }
  82.    };
  83.   };
  84.  
  85.  
  86.   {
  87.    int order = 102;
  88.    int uplo = 122;
  89.    int trans = 111;
  90.    int N = 1;
  91.    int K = 2;
  92.    float alpha = -0.3f;
  93.    float beta = 0.1f;
  94.    float A[] = { 0.655f, 0.45f };
  95.    int lda = 1;
  96.    float B[] = { 0.16f, -0.747f };
  97.    int ldb = 1;
  98.    float C[] = { 0.576f };
  99.    int ldc = 1;
  100.    float C_expected[] = { 0.19641f };
  101.    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  102.    {
  103.      int i;
  104.      for (i = 0; i < 1; i++) {
  105.        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1625)");
  106.      }
  107.    };
  108.   };
  109.  
  110.  
  111.   {
  112.    int order = 101;
  113.    int uplo = 121;
  114.    int trans = 112;
  115.    int N = 1;
  116.    int K = 2;
  117.    float alpha = -0.3f;
  118.    float beta = 0.0f;
  119.    float A[] = { 0.259f, -0.334f };
  120.    int lda = 1;
  121.    float B[] = { -0.911f, -0.426f };
  122.    int ldb = 1;
  123.    float C[] = { 0.432f };
  124.    int ldc = 1;
  125.    float C_expected[] = { 0.056199f };
  126.    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  127.    {
  128.      int i;
  129.      for (i = 0; i < 1; i++) {
  130.        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1626)");
  131.      }
  132.    };
  133.   };
  134.  
  135.  
  136.   {
  137.    int order = 101;
  138.    int uplo = 122;
  139.    int trans = 112;
  140.    int N = 1;
  141.    int K = 2;
  142.    float alpha = -0.3f;
  143.    float beta = 0.0f;
  144.    float A[] = { -0.765f, 0.7f };
  145.    int lda = 1;
  146.    float B[] = { 0.487f, 0.768f };
  147.    int ldb = 1;
  148.    float C[] = { 0.836f };
  149.    int ldc = 1;
  150.    float C_expected[] = { -0.099027f };
  151.    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  152.    {
  153.      int i;
  154.      for (i = 0; i < 1; i++) {
  155.        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1627)");
  156.      }
  157.    };
  158.   };
  159.  
  160.  
  161.   {
  162.    int order = 102;
  163.    int uplo = 121;
  164.    int trans = 112;
  165.    int N = 1;
  166.    int K = 2;
  167.    float alpha = -0.3f;
  168.    float beta = 0.0f;
  169.    float A[] = { -0.584f, 0.056f };
  170.    int lda = 2;
  171.    float B[] = { 0.928f, -0.101f };
  172.    int ldb = 2;
  173.    float C[] = { -0.529f };
  174.    int ldc = 1;
  175.    float C_expected[] = { 0.328565f };
  176.    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  177.    {
  178.      int i;
  179.      for (i = 0; i < 1; i++) {
  180.        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1628)");
  181.      }
  182.    };
  183.   };
  184.  
  185.  
  186.   {
  187.    int order = 102;
  188.    int uplo = 122;
  189.    int trans = 112;
  190.    int N = 1;
  191.    int K = 2;
  192.    float alpha = -0.3f;
  193.    float beta = 0.0f;
  194.    float A[] = { 0.25f, 0.8f };
  195.    int lda = 2;
  196.    float B[] = { 0.489f, -0.642f };
  197.    int ldb = 2;
  198.    float C[] = { 0.322f };
  199.    int ldc = 1;
  200.    float C_expected[] = { 0.23481f };
  201.    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  202.    {
  203.      int i;
  204.      for (i = 0; i < 1; i++) {
  205.        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1629)");
  206.      }
  207.    };
  208.   };
  209.  
  210.  
  211.   {
  212.    int order = 101;
  213.    int uplo = 121;
  214.    int trans = 111;
  215.    int N = 1;
  216.    int K = 2;
  217.    double alpha = 0;
  218.    double beta = 0;
  219.    double A[] = { 0.591, 0.21 };
  220.    int lda = 2;
  221.    double B[] = { -0.718, -0.579 };
  222.    int ldb = 2;
  223.    double C[] = { -0.856 };
  224.    int ldc = 1;
  225.    double C_expected[] = { -0.0 };
  226.    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  227.    {
  228.      int i;
  229.      for (i = 0; i < 1; i++) {
  230.        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1630)");
  231.      }
  232.    };
  233.   };
  234.  
  235.  
  236.   {
  237.    int order = 101;
  238.    int uplo = 122;
  239.    int trans = 111;
  240.    int N = 1;
  241.    int K = 2;
  242.    double alpha = 0;
  243.    double beta = 0;
  244.    double A[] = { -0.971, -0.824 };
  245.    int lda = 2;
  246.    double B[] = { -0.227, 0.457 };
  247.    int ldb = 2;
  248.    double C[] = { 0.521 };
  249.    int ldc = 1;
  250.    double C_expected[] = { 0.0 };
  251.    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  252.    {
  253.      int i;
  254.      for (i = 0; i < 1; i++) {
  255.        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1631)");
  256.      }
  257.    };
  258.   };
  259.  
  260.  
  261.   {
  262.    int order = 102;
  263.    int uplo = 121;
  264.    int trans = 111;
  265.    int N = 1;
  266.    int K = 2;
  267.    double alpha = 0;
  268.    double beta = 0;
  269.    double A[] = { -0.274, 0.583 };
  270.    int lda = 1;
  271.    double B[] = { 0.668, -0.83 };
  272.    int ldb = 1;
  273.    double C[] = { 0.907 };
  274.    int ldc = 1;
  275.    double C_expected[] = { 0.0 };
  276.    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  277.    {
  278.      int i;
  279.      for (i = 0; i < 1; i++) {
  280.        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1632)");
  281.      }
  282.    };
  283.   };
  284.  
  285.  
  286.   {
  287.    int order = 102;
  288.    int uplo = 122;
  289.    int trans = 111;
  290.    int N = 1;
  291.    int K = 2;
  292.    double alpha = 0;
  293.    double beta = 0;
  294.    double A[] = { -0.512, -0.436 };
  295.    int lda = 1;
  296.    double B[] = { -0.443, -0.259 };
  297.    int ldb = 1;
  298.    double C[] = { -0.667 };
  299.    int ldc = 1;
  300.    double C_expected[] = { 0.0 };
  301.    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  302.    {
  303.      int i;
  304.      for (i = 0; i < 1; i++) {
  305.        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1633)");
  306.      }
  307.    };
  308.   };
  309.  
  310.  
  311.   {
  312.    int order = 101;
  313.    int uplo = 121;
  314.    int trans = 112;
  315.    int N = 1;
  316.    int K = 2;
  317.    double alpha = -0.3;
  318.    double beta = 0;
  319.    double A[] = { 0.741, -0.341 };
  320.    int lda = 1;
  321.    double B[] = { 0.743, -0.315 };
  322.    int ldb = 1;
  323.    double C[] = { -0.776 };
  324.    int ldc = 1;
  325.    double C_expected[] = { -0.3947868 };
  326.    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  327.    {
  328.      int i;
  329.      for (i = 0; i < 1; i++) {
  330.        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1634)");
  331.      }
  332.    };
  333.   };
  334.  
  335.  
  336.   {
  337.    int order = 101;
  338.    int uplo = 122;
  339.    int trans = 112;
  340.    int N = 1;
  341.    int K = 2;
  342.    double alpha = -0.3;
  343.    double beta = 0;
  344.    double A[] = { 0.03, 0.175 };
  345.    int lda = 1;
  346.    double B[] = { -0.832, 0.291 };
  347.    int ldb = 1;
  348.    double C[] = { 0.281 };
  349.    int ldc = 1;
  350.    double C_expected[] = { -0.015579 };
  351.    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  352.    {
  353.      int i;
  354.      for (i = 0; i < 1; i++) {
  355.        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1635)");
  356.      }
  357.    };
  358.   };
  359.  
  360.  
  361.   {
  362.    int order = 102;
  363.    int uplo = 121;
  364.    int trans = 112;
  365.    int N = 1;
  366.    int K = 2;
  367.    double alpha = -0.3;
  368.    double beta = 0;
  369.    double A[] = { 0.476, 0.496 };
  370.    int lda = 2;
  371.    double B[] = { -0.626, -0.159 };
  372.    int ldb = 2;
  373.    double C[] = { -0.964 };
  374.    int ldc = 1;
  375.    double C_expected[] = { 0.226104 };
  376.    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  377.    {
  378.      int i;
  379.      for (i = 0; i < 1; i++) {
  380.        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1636)");
  381.      }
  382.    };
  383.   };
  384.  
  385.  
  386.   {
  387.    int order = 102;
  388.    int uplo = 122;
  389.    int trans = 112;
  390.    int N = 1;
  391.    int K = 2;
  392.    double alpha = -0.3;
  393.    double beta = 0;
  394.    double A[] = { -0.489, 0.611 };
  395.    int lda = 2;
  396.    double B[] = { -0.285, -0.673 };
  397.    int ldb = 2;
  398.    double C[] = { -0.11 };
  399.    int ldc = 1;
  400.    double C_expected[] = { 0.1631028 };
  401.    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  402.    {
  403.      int i;
  404.      for (i = 0; i < 1; i++) {
  405.        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1637)");
  406.      }
  407.    };
  408.   };
  409.  
  410.  
  411.   {
  412.    int order = 101;
  413.    int uplo = 121;
  414.    int trans = 111;
  415.    int N = 1;
  416.    int K = 2;
  417.    float alpha[2] = {1.0f, 0.0f};
  418.    float beta[2] = {0.0f, 0.1f};
  419.    float A[] = { 0.796f, 0.872f, -0.919f, 0.748f };
  420.    int lda = 2;
  421.    float B[] = { -0.945f, 0.915f, -0.252f, -0.276f };
  422.    int ldb = 2;
  423.    float C[] = { 0.07f, -0.957f };
  424.    int ldc = 1;
  425.    float C_expected[] = { -2.12843f, -0.054104f };
  426.    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  427.    {
  428.      int i;
  429.      for (i = 0; i < 1; i++) {
  430.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1638) real");
  431.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1638) imag");
  432.      };
  433.    };
  434.   };
  435.  
  436.  
  437.   {
  438.    int order = 101;
  439.    int uplo = 122;
  440.    int trans = 111;
  441.    int N = 1;
  442.    int K = 2;
  443.    float alpha[2] = {1.0f, 0.0f};
  444.    float beta[2] = {0.0f, 0.1f};
  445.    float A[] = { 0.984f, 0.526f, 0.284f, 0.806f };
  446.    int lda = 2;
  447.    float B[] = { -0.509f, -0.178f, 0.188f, -0.221f };
  448.    int ldb = 2;
  449.    float C[] = { -0.388f, 0.795f };
  450.    int ldc = 1;
  451.    float C_expected[] = { -0.43092f, -0.747044f };
  452.    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  453.    {
  454.      int i;
  455.      for (i = 0; i < 1; i++) {
  456.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1639) real");
  457.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1639) imag");
  458.      };
  459.    };
  460.   };
  461.  
  462.  
  463.   {
  464.    int order = 102;
  465.    int uplo = 121;
  466.    int trans = 111;
  467.    int N = 1;
  468.    int K = 2;
  469.    float alpha[2] = {1.0f, 0.0f};
  470.    float beta[2] = {0.0f, 0.1f};
  471.    float A[] = { -0.16f, 0.628f, -0.06f, -0.645f };
  472.    int lda = 1;
  473.    float B[] = { 0.846f, 0.545f, 0.032f, 0.493f };
  474.    int ldb = 1;
  475.    float C[] = { -0.041f, -0.621f };
  476.    int ldc = 1;
  477.    float C_expected[] = { -0.26101f, 0.783636f };
  478.    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  479.    {
  480.      int i;
  481.      for (i = 0; i < 1; i++) {
  482.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1640) real");
  483.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1640) imag");
  484.      };
  485.    };
  486.   };
  487.  
  488.  
  489.   {
  490.    int order = 102;
  491.    int uplo = 122;
  492.    int trans = 111;
  493.    int N = 1;
  494.    int K = 2;
  495.    float alpha[2] = {1.0f, 0.0f};
  496.    float beta[2] = {0.0f, 0.1f};
  497.    float A[] = { -0.478f, -0.556f, 0.519f, 0.177f };
  498.    int lda = 1;
  499.    float B[] = { -0.946f, 0.423f, -0.859f, 0.736f };
  500.    int ldb = 1;
  501.    float C[] = { -0.54f, -0.035f };
  502.    int ldc = 1;
  503.    float C_expected[] = { 0.226066f, 1.05345f };
  504.    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  505.    {
  506.      int i;
  507.      for (i = 0; i < 1; i++) {
  508.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1641) real");
  509.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1641) imag");
  510.      };
  511.    };
  512.   };
  513.  
  514.  
  515.   {
  516.    int order = 101;
  517.    int uplo = 121;
  518.    int trans = 112;
  519.    int N = 1;
  520.    int K = 2;
  521.    float alpha[2] = {0.0f, 0.0f};
  522.    float beta[2] = {1.0f, 0.0f};
  523.    float A[] = { -0.582f, 0.09f, -0.176f, 0.784f };
  524.    int lda = 1;
  525.    float B[] = { 0.687f, -0.859f, 0.945f, 0.756f };
  526.    int ldb = 1;
  527.    float C[] = { -0.663f, -0.186f };
  528.    int ldc = 1;
  529.    float C_expected[] = { -0.663f, -0.186f };
  530.    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  531.    {
  532.      int i;
  533.      for (i = 0; i < 1; i++) {
  534.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1642) real");
  535.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1642) imag");
  536.      };
  537.    };
  538.   };
  539.  
  540.  
  541.   {
  542.    int order = 101;
  543.    int uplo = 122;
  544.    int trans = 112;
  545.    int N = 1;
  546.    int K = 2;
  547.    float alpha[2] = {0.0f, 0.0f};
  548.    float beta[2] = {1.0f, 0.0f};
  549.    float A[] = { 0.231f, -0.452f, -0.112f, -0.837f };
  550.    int lda = 1;
  551.    float B[] = { -0.258f, 0.464f, -0.224f, 0.893f };
  552.    int ldb = 1;
  553.    float C[] = { -0.448f, 0.046f };
  554.    int ldc = 1;
  555.    float C_expected[] = { -0.448f, 0.046f };
  556.    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  557.    {
  558.      int i;
  559.      for (i = 0; i < 1; i++) {
  560.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1643) real");
  561.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1643) imag");
  562.      };
  563.    };
  564.   };
  565.  
  566.  
  567.   {
  568.    int order = 102;
  569.    int uplo = 121;
  570.    int trans = 112;
  571.    int N = 1;
  572.    int K = 2;
  573.    float alpha[2] = {0.0f, 0.0f};
  574.    float beta[2] = {1.0f, 0.0f};
  575.    float A[] = { 0.115f, 0.178f, -0.193f, -0.491f };
  576.    int lda = 2;
  577.    float B[] = { 0.545f, -0.665f, 0.979f, -0.4f };
  578.    int ldb = 2;
  579.    float C[] = { 0.522f, 0.712f };
  580.    int ldc = 1;
  581.    float C_expected[] = { 0.522f, 0.712f };
  582.    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  583.    {
  584.      int i;
  585.      for (i = 0; i < 1; i++) {
  586.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1644) real");
  587.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1644) imag");
  588.      };
  589.    };
  590.   };
  591.  
  592.  
  593.   {
  594.    int order = 102;
  595.    int uplo = 122;
  596.    int trans = 112;
  597.    int N = 1;
  598.    int K = 2;
  599.    float alpha[2] = {0.0f, 0.0f};
  600.    float beta[2] = {1.0f, 0.0f};
  601.    float A[] = { -0.725f, -0.808f, -0.244f, 0.145f };
  602.    int lda = 2;
  603.    float B[] = { 0.447f, -0.413f, -0.226f, -0.585f };
  604.    int ldb = 2;
  605.    float C[] = { -0.531f, 0.227f };
  606.    int ldc = 1;
  607.    float C_expected[] = { -0.531f, 0.227f };
  608.    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  609.    {
  610.      int i;
  611.      for (i = 0; i < 1; i++) {
  612.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1645) real");
  613.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1645) imag");
  614.      };
  615.    };
  616.   };
  617.  
  618.  
  619.   {
  620.    int order = 101;
  621.    int uplo = 121;
  622.    int trans = 111;
  623.    int N = 1;
  624.    int K = 2;
  625.    double alpha[2] = {-0.3, 0.1};
  626.    double beta[2] = {-0.3, 0.1};
  627.    double A[] = { 0.337, -0.737, -0.993, 0.69 };
  628.    int lda = 2;
  629.    double B[] = { -0.39, -0.836, -0.32, 0.368 };
  630.    int ldb = 2;
  631.    double C[] = { 0.844, -0.763 };
  632.    int ldc = 1;
  633.    double C_expected[] = { 0.3494384, 0.5248712 };
  634.    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  635.    {
  636.      int i;
  637.      for (i = 0; i < 1; i++) {
  638.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1646) real");
  639.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1646) imag");
  640.      };
  641.    };
  642.   };
  643.  
  644.  
  645.   {
  646.    int order = 101;
  647.    int uplo = 122;
  648.    int trans = 111;
  649.    int N = 1;
  650.    int K = 2;
  651.    double alpha[2] = {-0.3, 0.1};
  652.    double beta[2] = {-0.3, 0.1};
  653.    double A[] = { 0.386, -0.465, 0.719, -0.378 };
  654.    int lda = 2;
  655.    double B[] = { 0.099, -0.879, 0.864, 0.141 };
  656.    int ldb = 2;
  657.    double C[] = { -0.599, -0.47 };
  658.    int ldc = 1;
  659.    double C_expected[] = { 0.1664126, 0.5082238 };
  660.    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  661.    {
  662.      int i;
  663.      for (i = 0; i < 1; i++) {
  664.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1647) real");
  665.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1647) imag");
  666.      };
  667.    };
  668.   };
  669.  
  670.  
  671.   {
  672.    int order = 102;
  673.    int uplo = 121;
  674.    int trans = 111;
  675.    int N = 1;
  676.    int K = 2;
  677.    double alpha[2] = {-0.3, 0.1};
  678.    double beta[2] = {-0.3, 0.1};
  679.    double A[] = { 0.914, 0.128, -0.262, -0.26 };
  680.    int lda = 1;
  681.    double B[] = { 0.431, 0.276, 0.75, 0.904 };
  682.    int ldb = 1;
  683.    double C[] = { 0.287, 0.537 };
  684.    int ldc = 1;
  685.    double C_expected[] = { -0.3532044, 0.0216788 };
  686.    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  687.    {
  688.      int i;
  689.      for (i = 0; i < 1; i++) {
  690.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1648) real");
  691.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1648) imag");
  692.      };
  693.    };
  694.   };
  695.  
  696.  
  697.   {
  698.    int order = 102;
  699.    int uplo = 122;
  700.    int trans = 111;
  701.    int N = 1;
  702.    int K = 2;
  703.    double alpha[2] = {-0.3, 0.1};
  704.    double beta[2] = {-0.3, 0.1};
  705.    double A[] = { 0.618, 0.72, 0.392, -0.737 };
  706.    int lda = 1;
  707.    double B[] = { 0.783, 0.531, 0.375, 0.203 };
  708.    int ldb = 1;
  709.    double C[] = { 0.058, -0.116 };
  710.    int ldc = 1;
  711.    double C_expected[] = { -0.3837348, -0.2968344 };
  712.    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  713.    {
  714.      int i;
  715.      for (i = 0; i < 1; i++) {
  716.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1649) real");
  717.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1649) imag");
  718.      };
  719.    };
  720.   };
  721.  
  722.  
  723.   {
  724.    int order = 101;
  725.    int uplo = 121;
  726.    int trans = 112;
  727.    int N = 1;
  728.    int K = 2;
  729.    double alpha[2] = {0, 1};
  730.    double beta[2] = {1, 0};
  731.    double A[] = { -0.372, -0.735, -0.711, 0.051 };
  732.    int lda = 1;
  733.    double B[] = { 0.257, 0.097, 0.338, -0.484 };
  734.    int ldb = 1;
  735.    double C[] = { -0.142, -0.197 };
  736.    int ldc = 1;
  737.    double C_expected[] = { -0.414766, -0.676886 };
  738.    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  739.    {
  740.      int i;
  741.      for (i = 0; i < 1; i++) {
  742.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1650) real");
  743.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1650) imag");
  744.      };
  745.    };
  746.   };
  747.  
  748.  
  749.   {
  750.    int order = 101;
  751.    int uplo = 122;
  752.    int trans = 112;
  753.    int N = 1;
  754.    int K = 2;
  755.    double alpha[2] = {0, 1};
  756.    double beta[2] = {1, 0};
  757.    double A[] = { 0.1, -0.878, 0.28, -0.381 };
  758.    int lda = 1;
  759.    double B[] = { -0.208, 0.309, -0.276, 0.123 };
  760.    int ldb = 1;
  761.    double C[] = { 0.483, -0.541 };
  762.    int ldc = 1;
  763.    double C_expected[] = { -0.22324, -0.10083 };
  764.    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  765.    {
  766.      int i;
  767.      for (i = 0; i < 1; i++) {
  768.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1651) real");
  769.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1651) imag");
  770.      };
  771.    };
  772.   };
  773.  
  774.  
  775.   {
  776.    int order = 102;
  777.    int uplo = 121;
  778.    int trans = 112;
  779.    int N = 1;
  780.    int K = 2;
  781.    double alpha[2] = {0, 1};
  782.    double beta[2] = {1, 0};
  783.    double A[] = { -0.918, 0.515, -0.985, 0.067 };
  784.    int lda = 2;
  785.    double B[] = { -0.034, 0.408, 0.66, -0.945 };
  786.    int ldb = 2;
  787.    double C[] = { -0.063, -0.018 };
  788.    int ldc = 1;
  789.    double C_expected[] = { -1.228982, -1.549386 };
  790.    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  791.    {
  792.      int i;
  793.      for (i = 0; i < 1; i++) {
  794.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1652) real");
  795.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1652) imag");
  796.      };
  797.    };
  798.   };
  799.  
  800.  
  801.   {
  802.    int order = 102;
  803.    int uplo = 122;
  804.    int trans = 112;
  805.    int N = 1;
  806.    int K = 2;
  807.    double alpha[2] = {0, 1};
  808.    double beta[2] = {1, 0};
  809.    double A[] = { 0.443, -0.009, -0.245, -0.008 };
  810.    int lda = 2;
  811.    double B[] = { 0.495, -0.239, -0.973, -0.032 };
  812.    int ldb = 2;
  813.    double C[] = { -0.85, -0.799 };
  814.    int ldc = 1;
  815.    double C_expected[] = { -0.660584, 0.111526 };
  816.    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  817.    {
  818.      int i;
  819.      for (i = 0; i < 1; i++) {
  820.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1653) real");
  821.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1653) imag");
  822.      };
  823.    };
  824.   };
  825.  
  826.  
  827. }
  828.